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Abstract 

It is generally expected that adding light sterile species would increase the effective number of 
neutrinos, N e g. In this paper we discuss a scenario that N e g can actually decrease due to the 
neutrino oscillation effect if sterile neutrinos have self-interactions. We specifically focus on the eV 
mass range, as suggested by the neutrino anomalies. With large self-interactions, sterile neutrinos 
are not fully thermalized in the early Universe because of the suppressed effective mixing angle or 
matter effect. As the Universe cools down, flavor equilibrium between active and sterile species 
can be reached after big bang nucleosynthesis (BBN) epoch, but leading to a decrease of N e g. In 
such a scenario, we also show that the conflict with cosmological mass bounds on the additional 
sterile neutrinos can be relaxed further when more light species are introduced. To be consistent 
with the latest Planck results, at least 3 sterile species are needed. 
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I. INTRODUCTION 


Although most neutrino experiments can be well described by the standard three light 
species paradigm, there have been several anomalies that indicate a new sterile state [1-6]. 
The new state should have an eV scale mass and a mixing angle with sin 2 29o ~ 0.01 [7, 8], 
With such parameters, sterile neutrinos can be copiously produced from oscillation with each 
species increasing the effective number of neutrinos N e g almost by one unit. This would be 
in tension with the cosmological bounds from cosmic microwave background (CMB) data [9], 
iVeff < 3.91 and < 0.59eV. 

The conflict can be resolved if the sterile species is only partially thermalized with N e g <C 
1. Partial thermalization can be realized in particle physics models where there are secret 
self-interactions 1 in the sterile neutrino sector [10-14], These self-interactions can induce 
large matter potentials, effectively suppress the mixing angle, and block sterile neutrino’s 
production from oscillation efficiently [10, 11]. 

However, the situation changed recently when authors in [15] argued that even if sterile 
neutrinos’ production are blocked at BBN time, as the Universe cools down, the new self¬ 
interaction will eventually equilibrate sterile neutrinos with the active ones. Then, it can 
be easily shown that in 3+1 scenario, 1/4 of the cosmic background neutrinos would be the 
heavy sterile ones, which is still in conflict with the above cosmological bounds. 

The novelty in this paper is that, we numerically solve the quantum kinetic equations 
for neutrino mixing, show and confirm that flavor equilibrium is indeed reached after BBN, 
see Fig. 2. We also propose that the above mentioned conflict can be reconciled easily in 
an extension that more than one sterile states are introduced. Of the n introduced self¬ 
interacting sterile species, only one has eV scale mass and mixes with the active neutrinos. 
In the early Universe, all of them are out of equilibrium, but can approach flavor equilibrium 
with active neutrinos after BBN era. So the current relic neutrinos are composed of 3+n 
species with equal flavors. As we shall show that N e g in the late Universe can actually 
decrease and the cosmological bounds can be evaded. To be consistent with the latest 
Planck results, we found that n — 3 is the minimal number of the introduced sterile species. 

This paper is organized as follows. In section II, we briefly review the effective number 
of neutrinos, N e g. In section III, we discuss the scenario that how N e g in BBN epoch can be 

1 In this paper, we do not discuss the case that large lepton asymmetry exists in the active neutrino sector. 
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different from its value in CMB time when there are secret self-interactions among sterile 
neutrinos. In section IV, we show how the cosmological mass bounds on sterile neutrino can 
be relaxed if more than one light species are introduced. Finally, we give our conclusion. 


II. OVERVIEW OF iV eff 

In this section, we shall review N e g briefly and establish the related conventions and 
definitions. 

First, let us recall the thermal history of neutrinos in the standard model (SM) at the 
early Universe. When the temperature of the thermal bath, T 7 , is much higher than MeV, 
active neutrinos, v a {a = e, /i, r), are in thermal equilibrium with other SM particles through 
clectroweak processes and have the same temperature as T 7 . Around 2 MeV, active neu¬ 
trinos are decoupled because electroweak interaction is not strong enough to keep them in 
equilibrium. Later, electrons/positrons annihilate but heat only the photons. Using the 
conservation of entropy density, one can obtain the temperature ratio after e ± annihilation, 



( 2 . 1 ) 


Afterwards, T Va /T 1 is constant in the standard cosmology. 
Agff is defined by the energy density ratio, 



( 2 . 2 ) 


where Pr stands for the total energy density of radiations. In the standard model, we have 
3 species of active neutrinos, z/ e , and u T , so N eS = 3, or 3.046 precisely if instantaeneous 
neutrino decoupling is relaxed. If there is some new physics, it might modify T Va /T 1 and/or 
contribute to pr as extra radiation. For example, dark matter may affect T Ua [16, 17] or 


T 7 [18], and new particles can decay to or contribute as equivalent neutrinos [19-21], In this 


paper, we only consider the extended models with sterile neutrinos. 

In the following discussions, we shall use u a (a = e,p,r) to denote the 3 active neutrinos, 
u s for the sterile species, and up or v for all of them if not otherwise specified. And v t will be 
referred as the f-th mass eigenstate with a mass m Vi . For the parameter space we focus on, 
^ 1 , 2,3 are mainly mixed states of active neutrinos v a and zys (i > 3) are mainly mixed sterile 
species. We shall also neglect the masses of z/^ 2,3 when considering the mass constraints on 
sterile neutrinos. 
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FIG. 1. Thermal history of active/sterile neutrinos. When the temperature is high, v s s are not 
in thermal equilibrium with v a s because of the suppression from a large matter potential. As the 
Universe cools down, equilibrium between active and sterile neutrinos could be reached. 

In principle, N eS can be a function of time or photon temperature. We define SN eS as 
the deviation from the standard value, with explicit time/temperature dependence, 

^ LtL =l(^Pl I 3 + »-«]' < 2 - 3 ) 

Pi 8 V T i J 

or 

«.»(*) = £ - 3 . < 2 ' 4 ) 

P V 1 IP a / 

where T® a stands for the neutrino temperature in the standard cosmology without new 
physics, T® a = T 7 before annihilation and = (4/11 ) 1//3 T 7 afterwards, and /3 runs 
through all active/sterile neutrinos. 

We shall keep in mind that 3 in Eqs. 2.3 and 2.4 is actually 3.046 precisely. But this little 
difference would not affect our later discussions and we shall use 3 throughout the paper. 


III. 


JV e b ff bn VS ATf™ b 


In this section, we discuss how self-interaction can affect the thermalization of sterile 
neutrinos. The essential picture is described in Fig. 1 where sterile neutrinos are only 
partially thermalized at/before BBN time, but flavor equilibrium, p Vs = p Va , is reached at 
later time. To be as general as possible, we introduce n sterile species and do not discuss the 
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specific particle physics models, but emphasize that new interaction for the sterile neutrinos 
is required. 

It should be pointed out that in the present discussion for illustrating our points, the 
minimal setup we need is that one of the n sterile species is mixed with active neutrinos 
and has eV-scale mass, and the rest may be massless and have negligible mixing to affect 
current neutrino experimental results. Moreover, sterile neutrinos do not have to mix with 
each other since all sterile neutrinos can be in flavor equilibrium through new interactions. 
Therefore, let us just simply assume there is only mixing between a 4-th neutrino and the 
active species, namely the mixing matrix has the following form, 


( V e ^ 


( U el U e2 U e3 U e 4 • ■ ■ \ 


( V1 \ 



Ufti Ufj,2 3 Ufj ,4 ■ ■ ■ 


v 2 

V T 

= 

Url U T 2 U T 3 I/7-4 


V3 

V s 


U s 1 U s 2 Ug 3 ■ ■ ■ 


V4 

\ '■ ) 


\ ■■ 1 ) 


\ / 


( 3 . 1 ) 


Further assume only 1-4 mixing and CP conservation, then the complete 4x4 mixing matrix 
is 
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(3.2) 


where e.g. c t j = cos 0 VJ and ,s t j = sin 9 ^, etc. Complete investigations with the above 
multiple-flavor mixing are quite involved numerically when solving the full quantum kinetic 
equations (QKEs) [22-25] and we refer to Refs. [26-28] for multiple-flavor analysis with 
non-interacting sterile neutrinos. For simplicity and without loss of generality, we work with 
only two neutrino states, v e -v s mixing with mass difference 8m 2 and mixing angle 9 0 = 6 * 14 , 
and shall pay our attention to eV sterile neutrinos with the parameter space suggested by 
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neutrino anomalies, 5m 2 ~ leV 2 and sin 2 26 o ~ 0.01. We shall keep in mind that v s produced 
from and i+’s oscillation could be equally important since \U^\ ~ \U T ^\ ~ 0.46|t/ e 4 1. 

In most of the previous discussions in the literature, sterile neutrinos are assumed to have 
no interaction, so they can only be produced by oscillations from the active ones. And the 
production rate and total amount depend on the mass difference, mixing angle and lepton 
asymmetry [29-40]. 

In case of sterile species having a secret self-interaction, parametrized by G \ = fix/Mx 
(similar to the Fermi’s constant Gp in SM), 


G x VsFiV s VsTjVs, Fij are products of 7 ^, 75 ,..., 


(3.3) 


the production will also depend on the strength of Gx- Models that can give rise to the 
above types of self-interactions can be found in Refs. [10-14, 41, 42], To calculate the N^ n , 
we use the modified version of LASAGNA [10, 43, 44] with gx = 0.1. 

For 1+1 scenario, we can parametrize the system with 2x2 Hermitian density matrix in 
terms of Pauli matrices, 

p=\fo(Po + P-f), (3.4) 

where / 0 = 1/ [e E ^ T + l) and P = (P x , P y , P z ). And the densities of active and sterile 
neutrino are given by 


P a = P 0 + P z = 2 ^,P S = P 0 -P Z 
Jo 

The kinetic equations governing P^s time evolution are 



Pa 

Ps 


V x Py + r a 
— Vx Py + r s 


Je, 


2 j eq,a 


/o 

r -,feq,s 

fo 


~Pa 


~Ps 


Px = —VzPy — DP X , 

Py = V Z P X ~ -V x (P a ~ P s ) — DP y , 


(3.5) 

(3.6) 

(3.7) 

(3.8) 

(3.9) 


where V) and T, are the potentials and scattering kernels [ 10 ] , respectively, 


t 1 &m 2 . n T f 5m 2 14vr 2 G F 4 14vr 2 G x 
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FIG. 2. Evolution of 5N e g as temperature T 7 decreases. SN^g n depends on the self-interaction 
strength Gx■ Black curve shows the non-interacting case, Gx = 0. The self-interaction can 
suppress the production of sterile neutrino at high temperature, but lead to flavor equilibrium at 
later time. Increasing the strength of self-interacting would delay the equilibrium time. 

From the QKEs, we can recover neutrino oscillation in vacuum if we take the non-interacting 
limits, Gp —t 0, Gx —t 0, 

P 0 = 0, P = VxP, 

which describe the precesses of P around V. Non-zero Tj’s effects are to repopulate different 
momentum modes to reach thermal distribution and P-terms would damp and shrink P. 
As we can see in the above equations, the introduced new self-interaction leads to T s ^ 0 
and contributes to 14 and D in the above QKEs. 

The introduced self-interactions for sterile neutrinos have two effects. One is to block the 
thermalization at high temperature. The other effect is the collisions that lead a scattering- 
induced decoherent production at later time. In Fig. 2, we show the non-interacting case 
with a black curve and compare it with two interacting cases, Gx = lOGp and Gx = 10 4 Gp. 
If Gx is large, for instance Gx = 10 4 GV, v s can only be partially thermalized with T Ug < T Ua 
when T 7 > 2MeV. However, if Gx is not large enough, for instance Gx = 10GV, it will block 
the thermalization first at high temperature but enhance the production of sterile neutrino 
at a later time even before BBN epoch, T 7 ~ 5MeV. The G\ = 10 4 GV case, however, 
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has shown that the equilibrium time is later than the Gx = 10 Gf case, less than IMeV. 
Generally, increasing G\ would delay the equilibrium time. 

A simplified picture to understand these two effects is to use the effective mixing angle, 

sin 2 2 On 


sin 2 6 L ff = 


(cos 200- HrV'ei ( ) 2 + sm 2 2So’ 


(3.10) 


where E is the energy of oscillating neutrino, 


5m 2 

~2E 


008 26*0 is usually called the vacuum 


oscillation term [45], and V e g = Uj M — V[^ EW is the matter potential. Simple analysis would 
give 

_ 

(3.11) 


K SM ... 
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F ET‘ v„ 
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x et‘, 


which highly depend on the temperature. A familiar case, MSW effect [46-48], happens 
when V e R = cos 29 0 5m 2 /2E, leading to a maximal mixing angle 0 e g = 7 t/ 4 . However, when 
the matter potential is much larger than the vacuum term, the mixing angle is effectively 
suppressed, 


sin 2 26* e ff <C sin 2 26*o, when \V e g\ 


5m 2 


2 E 


COS 2 Or 


(3.12) 


which can efficiently block the production of sterile neutrinos. The value of T Va to block 
production at BBN time can be roughly estimated as follows: 


V e 


eff 


Gx t > Sm 2 


T 

-J- V 


M 2 X 


2 E 


MeV 


> 


5m 2 Ml \ 1/4 


2 E 2 G x 

Take Gx ~ lO 4 ^^ and M\ ~ 1 GeV, we get T Vs ~ 10 _3 MeV around BBN time. This 
is what wee see the smallness of T Us all the way to BBN time from the dot-dashed line in 
Fig. 2. 

As the Universe cools down, the matter potential V e g gets smaller very quickly. When 

Sm 2 

\Ves\ < yypr cos26*o , matter effect can be neglected and 6 c q ~ 6* 0 . Again, we discuss in the 
2 E 

simplified framework of two-flavor case, v a -v s or zq-z /2 in the mass eigenstates, 

= cos 6 * e ff ui — sin 6 * eff u 2 , v s = sin 6 * e fr + cos 6 * e ff z/ 2 - 


zy, 


Before BBN, because of the highly suppressed 6* e fr, there are mostly u a or u\ neutrinos in 
the thermal bath, and a small amount of u s or v 2 states. After 6 e g is not suppressed any 
more, zq has the scattering process to produce u 2 through zq + v 2 —>■ v 2 + v 2 with rate 
T ~ G\Tl sin 2 26q, and zq + ui —> u 2 + z / 2 with rate T ~ G\T^ a sin 4 2 0 o . If T is larger than 
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FIG. 3. vs 5N^ h . (Left panel)We choose several cases for the number of sterile species, 

as indicated by n. For each case, the solid curve shows while the dashed one gives 5N^°. 

Sizable differences can arise in the low T Ua /T® a region. (Right panel)<5Ag™ b as function of SN^q d . 
The solid black line is for SN^ h = 5N^g n in non-interacting case, and from up to down, dashed lines 
correspond to n = 1, 2, 3, 6. Dot-dashed orange lines mark the current bounds from Planck [49]. 


Hubble parameter H = 


87 tG 


p R , then Ui and u 2 , or u a and v si will reach equal numbers 


quickly. Since the sterile neutrinos’ self-interaction will induce rapidly elastic scattering, 
v s v s —s- u s u s , it can redistribute momenta among sterile neutrinos. As long as the scattering 
rate T x ~ G x T^ s is larger than the Hubble parameter H ( G is the Newton’s constant), sterile 
neutrinos will soon reach the Fermi-Dirac distribution, leading to the flavor equilibrium, 
n Us = n Ua - For G x = 10 8 Gi? and M x = 1.2MeV it was estimated that the equilibrium 
would be approached around T 7 ~ 40KeV [15], although numerically being a formidable 
challenge for such a large interaction. 


Before a detailed discussion on N e g at or after CMB time, we should note that the 
above investigation only took into account the sterile neutrinos produced from oscillations. 
However, oscillation is not the only contributing process. Generally, we expect there is a 
whole dark sector accompanying with self-interacting sterile neutrinos. Once the whole dark 
sector is decoupled from the SM thermal both, entropy is transfered to sterile neutrinos (and 
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other particles in thermal equilibrium with sterile neutrinos). We may call this part as the 
“primordial” portion (denoted as SN^g) which depends on the physical degrees of freedom in 
the dark sector and decoupling temperature, therefore model-dependent, see, for example, 
Refs. [12, 13] for concrete models. Since at high temperature oscillation is effectively blocked 
by the large matter potential, the effects of these primordial sterile neutrinos are to change 
the initial condition for QKEs and lift up the curves at high temperature by 5N® S in Fig. 2. 
Hence, if not stated explicitly, we shall treat T Us at BBN time as a free parameter in the 
rest of our discussion. 

Now, we are in a position to discuss the effect on N e $ at or after CMB time. Assume 
there are n sterile species with common temperature T Ug , if sterile and active neutrinos reach 
the equilibrium when they are still relativistic, they would have the same temperature T v 
determined by the conservation of neutrino number or entropy density, 

3 x (T°y + n x Tl = (3 + n) x (3.13) 


where n is the number of sterile species that have self-interactions 2 . 

With the new temperature T u , we can calculate SN e g at CMB time with Eq. 2.4, 

3-1 4/3 

3 + nx lp L l “3. (3-14) 


5N™ h = (3 + n)" 1/3 x 


in comparison with 

SNT = » x (f/) • (3-16) 

If sterile neutrinos were fully thermalized at BBN time, we would have T Ua = , and 

Eq. 3.14 gives the same result as Eq. 3.15 does. However, for partially thermalized u s , 
T Us < T° a . As shown in Fig. 3, it is evident that SN^ h < SN^ n , and that the difference 
can be significant in the low T Va /T® region and it increases as n gets bigger. An interesting 
observation is that SN^ h can even be negative for small values of T Us /T® a . If future experi¬ 
ment data indicate a deficit of N e g, it would be natural to consider the scenario that active 
neutrinos are mixing with self-interacting sterile species. 

In the right panel of Fig. 3, we show SN^ h as functions of ^A^ 11 for different n. The 
solid black line is for SN^ lh = (fA^ 11 in non-interacting case, and from up to down, dashed 


We note that n=l case has been discussed in the latest version of Ref. [15] whose results agree with ours. 
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lines correspond to n = 1, 2, 3, 6, respectively. Dot-dashed orange lines mark the current 
95% CL preferred ranges for N e g from Planck [49], 

N e ff = 3 .III 0 J 7 H e+Planck TT+lowP. (3.16) 

Since 5N b g n > 0 in our scenario, the available parameter space is inside the region with 
arrows indicated. 


IV. COSMOLOGICAL NEUTRINO MASS BOUNDS 


In this section, we show how the conflict between eV sterile neutrino and cosmological 
mass bounds can be relaxed when more than one light sterile species are introduced. 

Cosmological bounds on the neutrino masses from the combination of CMB, large scale 
structure and distance measurements are constraining the following effective quantity [9], 


eff _ n Vi m Vi 

"C - 0 


n 


Va 


( fit ) - 94.1eV x nji 2 , 


(4.1) 


where n Ui stand for number density, n b a for the value of active neutrino in standard 
cosmology, Pl v h 2 accounts for its energy density fraction in the Universe. After the later 
flavor equilibrium discussed above, all neutrinos share the same temperature, T Vi = T v . 
Using the same minimal setup as we did in last section, we assume only one sterile state has 
eV-scale mass and all the others are almost massless. So we can reduce the above summation 
only over the heaviest one, i = 4, 

T 

J- v 


eff 

— \ T 0 
Va 


m ,. 


(4.2) 


In Fig. 4, we plot how changes with T Vi /T b a in four cases, n — 1, 2, 3, 6. When more 
light states are added with but fixed 5N e g, the individual number density of each species is 
decreased. Therefore, the total number of the heaviest state u 4 is reduced and mf, s then gets 
smaller correspondingly. Other light states are just radiations and red-shifted, contributing 
only negligibly in late Universe. 

Now, we compare with the cosmological bounds. We should note that cosmological 
bounds on sterile neutrino mass and abundance depend on the cosmological models and the 
chosen data set [50-58] 3 . Varying N^ b and only, and using the Planck+WP+HighL 
3 Currently, there is a very loose constraint on self-interaction of active neutrinos from CMB data [58], 


around 


10 8 Gf. The bound on self-interaction of sterile neutrino can be inferred and should be similar. 
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FIG. 4. (Left panel)m® ff vs n. is a model-dependent quantity. We show how it changes as 
T Ua /T® a varies in four cases with different n. m Vi = leV is assumed to be dominant on the mass. 
(Right panel)m® ff as functions of 5N^° and n. Parameter space inside region marked with arrows 
is still allowed. 

data combination, the latest result from Planck collaboration was able to give bounds with 
95% CL [49], 

2.53 < iV e 7 b < 3.7, mf < 0.52eV. (4.3) 

As we show in right panel of Fig. 4, if we the face value of the above constraint 4 , then n — 3 
is on the intersect point with marginal status and is the minimal number of introduced 
sterile species. This amusing accidental agreement recovers the symmetry between active 
and sterile neutrinos. If future experimental pushed the upper limit further stringent, from 
the trend shown in the four cases of Fig. 4, it is easy to introduce more light sterile states in 
the discussed scenario to relax the cosmological bounds. When putting the lower bounds on 
TVeff, we should be aware of the assumption that no other relativistic particles contribute as 
radiations. In cases where there are quite a mount of massless particles such as Goldstone 
or Majaron particles, the lower bounds on IV e fr then do not apply and n < 3 will be allowed. 


4 Note that the above constraint is only intended for non-interacting sterile neutrinos. Self-interactions may 
change the “free-streaming” scale, see Ref. [59] for example. Since no analysis with self-interacting sterile 
neutrino is available, here we only took the face value from Planck [49]. 
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V. CONCLUSION 


In this paper, we have discussed a scenario that eV sterile neutrinos are partially ther- 
malized before BBN era but equilibrated with active ones in later time. A mechanism 
to realize such a scenario is to introduce secret self-interactions for sterile neutrinos. The 
self-interactions can induce large matter potentials at high temperature, suppress the mix¬ 
ing angle and block the production of sterile neutrinos from oscillations. They can also 
lead to a rapid scattering-induced decoherent production of sterile neutrinos at later times 
before CMB. When flavor equilibrium between active and sterile species is approached, it 
surprisingly leads to a decrease of N eS . 

We also discussed how the conflict with cosmological neutrino mass bounds can be re¬ 
laxed in this scenario. We have found that, the more light sterile species we add, the less 
constrained they would be. If we take the latest Planck bounds [49], < 3.7 and 

< 0.52eV, at least three sterile species are needed to evade such a constraint. 
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